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ABSTRACT 


This  report  describes  an  Integrated  package  for  evaluating  the  dynamic 
effects  of  driving  a vehicle  across  any  described  section  of  terrain. 

The  model  Incorporated  Includes  only  the  pitch  and  roll  motions  of  the 
chasls  (or  hull).  The  package  Is  composed  of  four  digital  computer 
programs.  Each  program  Is  designed  to  be  run  Interactively  and  does 
not  require  that  the  user  have  any  programming  experience.  The  programs 
provide  the  capability  to  create  vehicle  data  files,  perform  simulations, 
and  specify  the  form  and  types  of  the  output  desired.  A general  description 
of  the  model  structure,  Its  capabilities  and  its  limitations  Is  included 
along  with  user  Instructions  for  running  each  program  in  the  package. 
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INTRODUCTION  AND  OBJECTIVE 


The  US  Army  Tank-Automotive  Command  has  been  involved  in  the 
simulation  of  tracked  and  wheeled  vehicles  for  over  twenty  years. 

The  objective  of  one  type  of  simulation  often  used  in  the  evaluation 
of  vehicle  concepts  has  been  to  determine  the  ride  performance  of  a 
vehicle  as  It  traverses  a section  of  terrain  or  a discrete  obstacle. 

The  model  generally  used  for  these  simulations  represents  the 
vehicle  as  a sprung  mass  (the  chasis  or  hull)  and  a separate  mass 
for  each  wheel.  Springs  and  dampers  were  represented  by  piecewise 
linear  curves.  The  tires  were  represented  also  as  a vertical  spring 
and  damper  between  the  wheel  center  and  the  ground. 

This  type  of  model  was  originally  developed  and  Implemented  on 
a set  of  analog  computers  with  a digital  computer  providing  the  for- 
cing function  during  the  simulation.  This  required  specially  trained 
engineers  to  set  up  and  run  the  simulations,  and  required  extensive 
program  checkout  for  each  configuration  to  ensure  the  validity  of  model 
implementation.  As  a result,  vehicle  evaluations  that  were  run  on 
this  model  were  very  expensive.  Due  to  the  extensive  checkout  required 
for  each  implementation,  these  simulations  appeared  to  be  quite  reason- 
able and  when  field  test  data  were  available,  they  tended  to  reinforce 
this  judqement. 

In  recent  years,  it  has  often  become  necessary  to  obtain  results 
in  a matter  of  a few  days.  Using  the  analog  computers  available, 
each  concept  took  about  two  weeks  to  complete  with  very  little  effective 
overlap  between  concepts.  With  the  spread  of  the  use  of  high  speed 
digital  computers  in  the  Army  community,  it  became  feasible  to  do 
large  simulations  on  digital  computing  machines.  With  this  in  mind,  a 
digital  counterpart  of  the  analog  model  was  developed  under  contract 
by  Service  Bureau  Corporation  and  was  made  operational  at  TARADCOM 
atout  10  years  ago.  In  the  ensuring  years,  countless  vehicles  and  vehicle 
ct ncepts  have  been  simulated  and  evaluated  using  evolutionary  forms 
of  this  (or  very  similar)  digital  models. 

Several  factors  have  persistently  arisen  over  the  years  that  have 
lessened  the  overall  utility  of  this  type  of  digital  model.  One  of 
these  factors  is  the  converting  of  vehicle  data  into  a computer  use- 
able  format.  The  person  performing  such  a task  must  be  familiar  with 
both  the  vehicle  data  and  the  computer  program  requirements.  Further- 
more, he  must  not  take  his  task  casually,  since  mistakes  may  not  be 
readily  visible  in  the  results  but  can  cause  great  repercussions  If 
discovered  at  a later  time. 
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Another  problem  has  been  that  different  users  may  want  different 
types  of  output  from  the  model.  This  would  necessitate  program 
chanqes  which  would  slow  down  the  project  and  increase  the  possibility 
of  errors  in  the  model,  since  adequate  program  checkout  time  is  seldom 
available. 

A third  retarding  factor  is  that  the  engineer  responsible  for  the 
work  has  had  no  opportunity  to  gain  personal  confidence  In  the  model 
being  used.  If  the  person  making  the  actual  computer  runs  for  him 
made  any  data  translation  or  transcription  errors,  or  if  any  programming 
errors  were  made  resulting  in  unreasonable  output,  the  engineer  tended 
to  lose  confidence  in  the  model.  With  a couple  of  bad  experiences  under 
his  belt,  he  would  soon  be  very  hesitant  to  use  the  model. 

During  the  past  two  years  one  evolutionary  form  of  this  model  has 
been  resurrected.  A preprocessor  and  a post-processor  program  have 
been  developed  which  when  used  in  conjunction  with  the  digital  sim- 
ulation model,  provide  a total  vehicle  ride  evaluation  package  which 
hopefully  will  eliminate  many  of  the  most  severe  problems  encountered 
in  the  past.  This  package  has  been  implemented  on  the  CDC  6500/66(10 
computer  complex  located  at  Picatinny  Arsenal  in  Dover,  New  Jersey.  The 
underlying  theme  used  in  the  development  of  the  packaqe  is  that  the 
ennineer  responsible  for  the  results  should  become  personally  involved 
with  the  model.  In  order  to  accomplish  this  objective,  the  package  was 
designed  to  be  run  interactively  from  a time-sharinq  terminal.  The 
vehicle  data  are  entered  interactively  in  ennineering-type  terms  in 
response  to  specific  requests.  Each  simulation  over  a particular  ter- 
rain at  a specified  velocity  is  also  run  interactively,  and  the  standard 
measures  of  driver  ride  comfort  are  obtained.  The  post  processor 
portion  of  the  package  can  be  run  interactively  as  well,  giving  the 
user  the  option  to  obtain  various  tabular  or  graphical  results  from 
any  signa1  generated  during  the  simulation  portion  of  the  package.  The 
full  capability  of  the  post  processor  is  available  only  from  a graphics 
terminal  having  vector  plotting  capability.  (A  Tektronix  4014  is  used.) 

The  remainder  of  this  report  is  concerned  with  a user  oriented 
description  of  this  packaqe.  Included  are  a description  of  the  model, 
some  model  verification  results,  and  a discussion  of  how  to  use  each 
component  of  the  packaqe.  Finally,  included  as  Appendix  I,  is  a 
rather  complete  example  of  what  can  be  accomplished  with  this  model. 

The  example  begins  with  the  vehicle  data  in  the  form  an  engineer  might 
receive  them  and  proceeds  through  the  sequence  of  programs  to  create 
a vehicle  data  file,  run  the  simulation,  and  ultimately  obtain  additional 
output  via  the  post  processor  portion  of  the  package. 
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f'ODFL  DESCRIPTION  AND  CAPABILITIES 


The  model  implemented  in  this  package  is  represented  pictor- 
ially  in  Figure  1.  It  represents  a wheeled  or  tracked  vehicle 
(Note  that  a tracked  vehicle  is  handled  as  a wheeled  vehicle  for  the 
purpose  of  ride  simulation).  The  pitch  and  vertical  motions  of  the 
sprunn  mass  (chassis  or  hull)  are  modeled  as  well  as  the  vertical 
motion  of  each  wheel.  The  standard  output  of  the  program  is  the  peak 
and  rms  vertical  accelerations  at  the  driver's  station.  It  also  cal- 
culates the  absorbed  power  of  the  driver  from  the  acceleration  sig- 
nal at  the  driver's  station.  Absorbed  power  is  a frequency  weighted 
rms  acceleration.  The  weighting  is  done  based  on  human  ride  tolerance 
to  the  various  frequencies,  and  is  explained  in  reference  1. 

The  data  necessary  to  run  this  program  are  the  wheel  and  CG 
locations,  vehicle  sprung  and  unsprung  mass,  sprung  pitch  moment 
of  inertia,  suspension  spring  and  damping  curves,  tire  spring 
and  damping  characteristics,  and  a description  of  the  terrain 
over  which  the  vehicle  is  to  be  simulated.  The  program  allows  for 
independently  sprung  wheels  or  for  a pair  of  wheels  to  be  joined 
in  a tandem  arrangement.  The  program  also  allows  the  input  of  the 
vertical  clearance  and  horizontal  location  of  one  possible  impact 
point  (e.o.,  a sprocket  or  bumper),  and  allows  the  inclusion  of  a fifth 
wheel  type  trailer  hitch  with  full  pitch  freedom. 

A typical  output  from  this  model  will  contain  all  the  data  shown 
in  Figure  2.  Additionally  a binary  output  file  is  created  which  con- 
tains all  the  relevant  displacements,  velocities,  and  accelerations 
required  to  reconstruct  the  entire  simulation.  The  various  signals 
recorded  in  this  file  may  be  further  analyzed  using  a separate  post 
processor  program.  A detailed  description  of  how  to  use  this  capability 
is  contained  in  the  User  Instruction  Section  under  the  program  headinq, 
BANALYZE. 


The  simulation  program  contains  a set  of  non-linear  differential 
equations.  These  equations  describe  the  accelerations  of  the  various 
bodies  being  modeled  and  are  as  follows: 

The  vertical  acceleration  of  the  sprung  center  of  mass  is: 
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where 


■ the  spring  rate  at  wheel  1 
Mq  • 1/2  the  sprung  mass  of  the  vehicle 
Yj  - the  hull  reference  position  of  wheel  1 
YWj  ■ vertical  position  of  wheel  1 
n ■ number  of  wheels  on  one  side  of  vehicle 
Dj  - the  damping  rate  at  wheel  1 

2 

g ■ the  gravitational  constant  • 32.16  ft/sec 

The  pitch  acceleration  of  the  hull  about  the  sprung  center  of  mass 
Is: 


where 


Lj  ■ the  horizontal  distance  from  the  CG  to  wheel  1 
J0  • 1/2  the  sprung  pitch  moment  of  Inertia  of  the  vehicle 
The  vertical  acceleration  of  wheel  1 (for  1 ■ 1 to  n): 


+ 
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where 


■ the  mass  of  wheel  1 
» the  terrain  elevation  at  wheel  1 

KW1  * the  spring  rate  of  wheel  1 

Dw  ■ the  damping  rate  of  wheel  1 

W1 

In  addition  to  the  above  differential  equations,  a pair  of  equations 
for  each  wheel  are  used  to  describe  the  hull  reference  point  position 
and  velocity  associated  with  each  wheel.  These  give  the  following 
equations  for  wheel  1: 

Yj  = Y„  •+  Li  si  r>  © 

Y-.  * f,  + 6 Li  Cos  © 

where 

Y0  ■ the  vertical  position  of’ the  CG 

6 * the  pitch  angle  of  the  vehicle  (front  end  up  Is  positive) 

The  simulation  Is  performed  by  the  simultaneous  Integration  of 
the  differential  equations  concurrently  with  the  solution  of  the 
necessary  auxiliary  equations.  This  numerical  Integration  Is  per- 
formed using  a variable  step  size  numerical  Integration  alogrlthm. 

A fourth  order  Adams -Moulton  predictor  corrector  alogrlthm  Is  em- 
ployed whenever  sufficient  past  values  are  available.  A fourth 
order  Runge-Kutta  algorithm  Is  used  for  starting  the  simulation  and 
is  also  used  for  restarting  If  step  size  changes  have  made  the 
appropriate  past  values  unavailable. 

Acceleration  values  are  saved  at  equal  1 ttervals  of  time  for 
subsequent  analysis.  Upon  completion  of  tl  simulation  portion 
of  the  program,  the  program  determines  the  peak  and  RMS  accelerations 
for  the  run,  measured  at  the  driver's  station  (dr  at  some  other 
designated  point  of  Interest). 

The  power  spectral  density  (PSD)  of  the  driver's  vertical  acc- 
eleration Is  calculated  and  passed  through  the  absorbed  power 


weighting  function,  and  Integrated.  The  resulting  value  is  the 
driver's  vertical  absorbed  power  In  watts.  This  value  Is  printed 
along  with  the  RMS  and  peak  accelerations. 

The  forcing  function  for  the  simulation  is  the  ground  profile 
that  would  be  driven  over  by  one  side  of  the  vehicle.  This  profile 
may  be  a surveyed  course,  a discrete  rigid  obstacle,  or  a mathematically 
generated  terrain  possessing  certain  statistical  characteristics. 

The  terrain  profile,  regardless  of  its  origin,  it  represented  by 
a series  of  elevations  at  known  horizonal  positions  along  the  course. 
These  known  profile  points  are  assumed  to  be  connected  by  straight 
line  segments.  This  profile  representation  is  modified  by  slowly 
rolling  a rigid  wheel  across  it  and  determining  the  lowest  path  that 
this  rigid  wheel  could  follow.  This  minimum  path  represents  the 
original  path,  modified  to  reflect  the  effects  of  tire  bridging.  When 
simulating  tracked  vehicles,  the  roadwheel  radius  plus  the  track  thick- 
ness is  used  as  the  effective  radius  of  the  rigid  wheel,  whereas  for 
wheeled  vehicles  the  rolling  radius  of  the  pneumatic  tire  Is  used. 

This  terrain  modification  is  performed  before  the  simulation  is  actually 
begun.  The  resulting  minimum  profile  is  used  as  the  forcing  function 
inDut  to  the  bottom  of  the  wheel  during  the  simulation.  It  should  be 
noted  that  this  terrain  modification  can  be  suppressed,  if  it  Is  so 
desired,  by  setting  the  tire  diameter  equal  to  zero  in  the  vehicle 
data  file. 

The  model  implementation  allows  a maximum  of  ten  wheels  per  side 
of  the  vehicle.  Each  wheel  must  have  a spring  associated  with  it 
arid  may  also  have  a shock  absorber.  In  addition,  any  pair  of 
adjacent  wheels  may  be  joined  together  in  a tandem  arrangement.  The 
tandem  is  free  to  rotate  about  a point  on  the  sprung  mass  and  is  re- 
presented by  its  pitch  moment  of  inertia,  Its  attachment  point  to  the 
sprung  mass,  and  the  span  between  the  two  wheels  on  the  tandem. 

The  spring  and  damping  curves  describing  the  suspension  system 
elements  are  represented  as  piecewise  linear  functions  with  a maximum 
of  nine  linear  segments  each.  The  deflection  axis  for  the  spring 
curves  is  arbitrary,  but  it  must  be  consistent  from  wheel  to  wheel. 

In  other  words,  the  zero  deflection  point  may  correspond  to  zero 
force,  or  to  the  static  force,  or  some  other  force.  The  effect  of 
being  consistent  in  the  deflection  axis  definition  for  all  wheels  can 
be  described  by  the  following  situation.  If  the  chasis  were  constrained 


to  remain  level,  and  each  wheel  were  positioned  with  the  same  value  of 
deflection  (regardless  of  the  value  chosen)  across  its  respective  spring, 
then  the  wheel  centerlines  would  all  lie  on  a single  horizonal  line. 

It  should  then  be  noted  that  if  all  wheels  have  zero  force  at  zero  deflection 
or  if  all  wheels  have  static  force  at  zero  deflection,  then  consistency 
is  guaranteed.  One  of  these  two  definitions  is  generally  used. 

An  example  of  an  inconsistent  definition  would  be  to  select  the 
jounce  (compression)  stop  point  as  the  zero  deflection  point.  If  this 
were  done  for  a vehicle  that  allowed  one  wheel  to  travel  nearer  to  the 
plane  formed  by  the  chasis  pitch  and  roll  axes  than  some  other  wheel, 
the  resulting  representation  would  not  be  as  intended. 

The  wheels  are  not  assumed  to  be  rigid  during  the  simulation.  In 
fact,  the  wheel  (or  tire)  vertical  spring  and  damping  characteristics 
along  with  the  unsprung  mass  of  each  wheel  are  required  to  simulated 
the  interaction  between  the  wheel  and  the  ground. 

The  vertical  spring  and  damping  characteristics  of  the  wheels  may 
differ  from  wheel  to  wheel.  These  characteristics  are  considered  to 
be  piecewise  linear  functions  and  may  have  a maximum  of  six  linear 
segments  to  describe  each  curve.  For  the  wheel  spring  curves,  the 
point  of  zero  deflection  must  coincide  with  the  point  of  zero  force. 

Negative  forces  cannot  be  generated  by  the  wheel  characteristics 
within  the  model.  In  other  words,  the  wheel  cannot  be  pulled  down  by 
the  profile,  but  it  can,  of  course,  leave  the  ground. 


Figure  2 Typical  Output  from  3VRj.De. 


MODEL  VALIDATION 


The  user  of  any  model  would  like  to  have  some  Idea  of  the  model's 
prediction  accuracy.  In  an  attempt  to  provide  some  degree  of  con- 
fidence in  this  model,  useahle  field  test  data  were  sought. 

One  such  set  of  results  were  obtained  from  the  M60A1  STB  (Standard 
Torsion  Bar)  tests  which  were  run  by  the  Waterways  Experiment  Station 
with  the  assistance  of  TARADCOM  personnel.  There  were  two  phases  of 
the  test.  The  first  phase  was  run  at  Fort  Knox  in  May  1976  when 
several  modifications  of  the  M60  were  tested.  The  second  phase  in- 
cluded the  MICV  as  well  as  the  M60  and  was  conducted  at  Aberdeen  Proving 
Ground  in  November  1976.  The  vehicles  were  run  over  terrains  of  varying 
roughness  and  the  driver's  absorbed  power  was  recorded  for  each  vehicle 
velocity. 

Figures  3 and  4 show  the  velocities  at  which  the  driver  experienced 
six  watts  absorbed  power  in  the  tests.  The  driver  comfort  limiting 
velocity  Is  plotted  as  a function  of  terrain  roughness  (quantified 
by  Inches  rtns)  for  the  M60A1  STB  and  for  the  MICV.  These  two  vehicles 
were  also  simulated  over  a variety  of  terrains  and  the  velocity  at  which 
the  simulated  driver  experienced  six  watts  of  absorbed  power  was  recorded. 
The  results  obtained  from  these  simulations  for  the  M60A1  STB  and  for 
the  MICV  are  also  plotted  In  Figures  3 and  4 for  comparison  with  the 
field  test  results. 
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USER  INSTRUCTIONS 


This  section  contains  Instructions  for  accessing  and  running  the 
various  programs  that  make  up  this  vehicle  dynamics  package.  The 
various  options  and  limitations  associated  with  each  program  are  discussed. 

These  programs  were  designed  to  be  run  from  a time-sharing  terminal 
and  should  not  require  any  pronramming  expertise  of  the  user.  The  user 
will,  of  course,  have  to  be  able  to  use  a time-sharing  terminal.  If  the 
user  can  gain  access  to  the  ARRADCOM  computer  system  through  a time- 
sharing terminal,  this  section  will,  hopefully,  provide  him  with  the 
instructions  necessary  to  run  the  various  programs  In  this  package. 


BCRET 


This  program  Was  written  to  allow  the  user  to  easily  create  a 
formatted  vehicle  data  tile  which  can  then  be  used  as  input  by  the 
simulation  portion  of  this  package.  The  program  is  designed  to  be 
run  interactively  from  a time-sharing  terminal.  After  the  user 
has  initiated  the  program,  various  requests  are  made  of  the  user  in 
terms  he  can  readily  understand. 

The  first  thing  the  user  is  asked  is  what  he  would  like  to  do. 

This  question  will  also  be  asked  each  time  a phase  of  the  program 
has  been  completed.  The  user  has  several  options  at  this  point.  He 
might  select  to  simply  stop  the  proqram  or  he  could  select  the  data 
entry  option.  If  this  latter  selection  is  made,  the  user  will  be 
asked  for  vehicle  data.  It  has  been  attempted  to  make  this  phase  of 
the  proqram  seem  as  natural  as  possible  to  the  engineer.  The  vehicle 
data  are  requested  in  the  form  the  enqineer  would  most  likely  have 
received  them.  The  entered  data  are  converted  to  the  form  and  format 
required  in  the  vehicle  data  file  and  the  file  is  written  on  Tape  8. 

The  user  may  also  select  the  data  check  option.  This  option  can 
only  be  employed  from  a graphics  terminal  and  Is  selected  by  typing 
an  8 or  9.  The  number  typed  corresponds  to  the  tape  number  from  which 
vehicle  data  are  to  be  read.  Tape  8 and  Tape  9 are  used  for  two 
different  forms  of  the  vehicle  data  either  of  which  can  be  used  with 
the  simulation  program.  Tape  8 corresponds  to  the  first  data  file 
name  in  the  program  run  statement.  The  output  of  the  data  entry  portion 
of  the  proqram  is  written  on  Tape  8. 

If  the  data  check  option  has  been  selected  by  typing  an  8,  the 
vehicle  represented  by  these  data  is  simulated  for  3.8  seconds  at 
5 mph  over  perfectly  flat  ground.  The  pitch  and  vertical  displacement 
of  the  sprung  mass  (chasis  or  hull)  center  of  gravity  are  plotted  as  a 
function  of  time.  These  plots  should  be  relatively  flat  if  the  data 
Is  reasonable  (note  the  converse  is  not  promised).  If  the  vehicle 
stabilizes  with  an  excessive  pitch  angle,  the  relative  travel  in  the 
various  suspension  components  will  be  altered.  A simple  example  might 
illustrate  this  problem  best.  Suppose  data  for  a vehicle  were  entered 
that  had  12  inches  of  wheel  travel  from  the  static  position  to  the 
jounce  stop  (jounce  represents  compression)  for  each  wheel.  Suppose 
also  that  the  vehicle's  sprung  center  of  gravity  was  behind  the 
geometric  center  of  the  wheels,  far  enough  to  cause  a static  pitch  angle 
of  0.05  radians  (about  2.85  degrees).  Assume  also  that  the  front  wheel 


was  8 feet  forward  of  the  vehicle's  sDrunq  CG.  Then  the  horizonal  dis- 
tance of  8 feet  swunq  through  the  pitch  anqle  of  0.05  radians  would 
cause  the  sprung  mass  at  the  front  wheel  station  to  rise  0.05  radians 
times  8 feet,  or  0.4  feet.  This  would  mean  that  the  front  wheel  would 
have  to  he  driven  up,  not  only  1 foot,  but  1.4  feet  before  the  bump 
stop  would  he  encountered.  Note  that  the  other  wheels  would  also  be 
affected  but  to  a lesser  extent.  In  this  example  the  wheels  forward  of 
the  CG  would  gain  jounce  travel  and  the  wheels  aft  of  the  CG  would  have 
their  jounce  travel  reduced. 

A technique -for  overcoming  this  difficulty  is  included  in  this  program. 
In  order  to  implement  this  technique  the  user  must  select  the  zero 
pitch  option.  The  data  entry  phase  must  have  been  completed  prior  to 
initiation  of  the  zero  pitch  phase.  In  order  to  level  the  vehicle  a 
simulation  is  performed  (again  for  3.8  seconds  at  5 mph),  over  level 
ground  using  the  vehicle  data  obtained  from  Tape  8.  Upon  completion  of 
this  simulation  the  vehicle  Is  assumed  to  have  stabilized.  The  force 
required  by  each  sprlnq  to  support  the  vehicle  in  this  stabilized  position 
is  recorded.  The  user  is  asked  how  much  static  to  jounce  stop  travel 
should  be  available  at  each  wheel  station.  The  spring  curve  associated  with 
each  wheel  is  then  translated,  first  along  the  deflection  axis  and  then 
along  the  force  axis.  The  resulting  curve  for  each  wheel  will  have  the 
appropriate  force  at  the  zero  deflection  point  to  support  the  vehicle 
in  a level  position  on  level  ground.  The  vehicle  data  with  these  modified 
spring  curves  are  then  written  on  Tape  9.  A flag  with  a value  of  1 Is 
written  in  column  60  of  card  2 to  indicate  that  the  vehicle  data  in  this 
file  have  been  generated  by  this  leveling  process.  The  simulation  program 
tests  this  flan  and  uses  a simpler  initialization  procedure  when  the  flag 
is  set  to  a value  of  one. 

The  zero  pitch  option  does  not  produce  any  plotted  output  and  therefore 
need  not  be  run  from  a Graphics  terminal.  The  only  immediate  indication 
to  the  user  that  the  leveling  process  has  been  completed  is  that  the 
option  request  will  be  made  once  again. 

After  leveling  the  vehicle  the  user  can  (if  he  is  at  a graphics 
terminal)  perform  a data  check  on  the  leveled  vehicle  data  file  in  Tape 
9 by  responding  to  the  option  request  with  a °.  This  data  check  is 
similar  to  the  one  described  previously  to  check  the  data  on  Tape  8. 

The  only  difference  is  in  the  method  of  initializing  the  vehicle. 

Certain  points  should  be  emphasized.  Tape  8 and  Tape  9 correspond 
to  the  two  file  names  associated  with  the  program  at  run  time  with  Tape  8 
beinn  the  first  and  Tape  9 the  second.  If  no  data  file  names  are  given, 
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then  the  local  file  name  for  Tape  ft  will  default  to  VEHIN  and  Tape  9 will 
default  to  VEHLEV.  If  BCRET  is  run  with  a leveled  vehicle  data  file 
created  at  some  earlier  time,  and  qiven  as  the  first  file  name,  then  a 
data  check  option  of  8 would  be  required  to  check  that  file. 

A couple  of  restrictions  should  be  reemphasized.  The  maximum  number 
of  senments  allowed  for  any  sprinq  or  dampinq  curve  is  9.  The  curves 
are  entered  as  an  initial  and  final  slope  and  breakooints  to  define  in- 
terior linear  segments.  Also,  the  number  of  wheels  on  one  side  of  the 
vehicle  cannot  exceed  10. 

Appendix  I contains  an  interactive  example  using  BCRET  to  enter 
the  vehicle  data.  A description  of  the  contents  of  the  vehicle  data 
file  is  qiven  in  Appendix  II. 

To  actually  run  this  program,  the  user  must  first  access  the  ARRADCOM 
computer  system  from  a time-sharing  terminal.  One  logged  into  the  sys- 
tem, the  user  must  type  the  following  two  commands  to  run  BCRET: 

ATTACH,  BCRET,  ID  - DRDTAZSA 

BCRET,  VEHIN,  VEHLEV 

In  the  latter  of  these  commands  the  two  parameters  VEHIN  and 
VEHLEV  may  be  chanqed  or  left  off  altogether.  If  omitted  the  default 
names  are  those  given  in  the  statement  above  (i.e.  VEHIN  and  VEHLEV). 

In  either  case  the  first  parameter  name  will  be  the  local  file  name 
containing  the  data  written  on  Tape  8 by  the  data  entry  option  of  the 
program.  The  second  parameter  name  will  be  the  local  file  name 
corresponding  to  Tape  9 (the  leveled  vehicle  data  file).  If  the  zero 
pitch  option  (to  level  the  vehicle)  is  not  selected,  then  nothing  will 
be  written  on  Tape  9. 
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BVRIDE 


This  program  will  perform  a dynamic  simulation  of  a vehicle  as  it 
traverses  a specified  terrain  profile.  Only  the  pitch  and  vertical 
degrees  of  freedom  are  modeled  for  the  sprung  mass  (hull  or  chasis)  and 
the  vertical  response  for  each  wheel.  The  model  as  implemented  in  this 
program  has  the  full  capability  described  In  the  section  of  this  report 
titled  Model  Description  and  Capabilities. 

This  simulation  can  be  run  from  any  time  sharing  terminal.  There 
is  no  graphical  output  generated  by  BVRIDE  and  only  the  vehicle  velocity 
is  entered  Interactively.  Before  running  the  simulation  the  user  must 
extend  the  time  limit  allotted  for  each  program  run  on  the  time  sharing 
system.  He  may  also  desire  to  inform  the  system  of  any  Increased  line 
width  his  terminal  may  have  {the  defaults  for  these  two  Items  are  7 
seconds  and  72  characters).  These  two  items  can  be  accomplished  when 
the  user  first  signs  on  the  system  and  will  remain  in  effect  until  the 
terminal  is  logged  out.  The  following  two  comnands  perform  these  functions 

FTL.500 

SCREEN , 1 32 

The  number  500  in  the  first  command  is  the  maximum  run  time  allotted. 

It  is  the  number  of  seconds  written  in  octal  (base  8)  notation. 

The  vehicle  data  must  be  supplied  in  a file  according  to  the  format 
specifications  spelled  out  in  Appendix  II.  This  file  can  be  created  by 
the  user  or  can  more  efficiently  be  generated  by  the  program  BCRET  which 
is  also  described  in  this  report.  The  only  feature  not  available  from 
BCRFT  is  the  allowance  for  non-linear  tire  characteristics.  If  the  user's 
data  requires  this  feature  the  proper  data  must  be  inserted  into  the 
vehicle  data  file  according  to  the  format  specifications  in  Appendix  II. 

The  terrain  profile  data  must  also  be  supplied  in  a separate  file. 

The  terrain  data  in  the  file  will  be  extended  sufficiently  in  both  dir- 
ections to  allow  the  vehicle  to  sit  on  level  ground  prior  to  crossing  the 
given  profile  and  to  return  to  level  ground  upon  completion  of  the  profile 
traverse. 

The  user  can  create  his  own  terrain  data  files  according  to  the  re- 
quired format  shown  In  Appendix  III.  The  other  alternative  is  to  use  one 
of  the  standard  terrains  already  on  the  system.  A catalog  of  these 


standard  terrains  is  incorporated  as  Appendix  IV.  Any  of  these  terrains 
may  be  accessed  by  selecting  the  appropriate  file  name  and  typing  the 
following  command: 

ATTACH,  filename,  ID=DRDTAZSA 

Once  the  user  has  a vehicle  data  file  and  a terrain  data  file  re- 
siding as  local  files,  he  is  ready  to  run  the  simulation.  The  two  com- 
mands necessary  to  run  this  program  are: 

ATTACH ,BVRI DE , ID- DRDTAZSA 

BVRIDE.VEH, TER, BINARY 

The  parameters  in  the  second  statement  represent  the  local  file 
names  corresponding  to  the  vehicle  data  file,  the  terrain  data  file, 
and  the  binary  output  file  respectively.  These  names  should  be  replaced 
with  the  local  file  names  the  user  has  employed. 

The  program  requests  only  the  vehicle  velocity  (in  MPH)  and  then 
proceeds  with  the  simulation  and  the  output.  The  printed  output  in- 
cludes a short  description  of  the  vehicle  and  terrain  used  for  the 
run.  The  vehicle  velocity  and  length  of  run  are  also  printed  out 
followed  by  the  simulation  results.  The  only  results  printed  are  the 
RMS  (root  means  square)  and  peak  vertical  acceleration  at  the  driver's 
station  followed  by  the  driver's  vertical  absorbed  power. 

When  the  simulation  is  complete,  the  user  may  run  the  post  processor 
proqram  RANALY7E  (also  described  within  this  report)  with  the  binary 
output  file  generated  by  BVRIDE.  With  this  technique  it  is  possible 
to  obtain  a wide  variety  of  types  of  output  from  this  simulation  package. 
If  one  chooses  to  run  the  post  processor  at  some  later  time,  the  user 
may  wish  to  make  use  of  the  system  CATALOG  command  to  save  the  binary 
output  file.  If  this  is  not  done,  and  another  simulation  is  run  using 
the  same  local  file  name  for  the  binary  output  file,  the  original  contents 
are  over  written  and  only  the  output  of  the  most  recent  simulation  will 
remain . 
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BANALYZE 


In  the  past  it  has  often  occured  that  different  users  of  the 
vehicle  dynamics  models  would  have  varying  objectives.  One  user  might 
be  quite  satisfied  with  simply  the  ride  comfort  level  as  a function  of 
vehicle  velocity  and  terrain  profile.  The  next  might  be  more  concerned 
with  the  force  levels  encountered  in  a spring  or  shock  absorber  duripg  a 
particular  run.  Changing  the  output  specifications  would  necessitate 
changes  to  the  programs  being  used.  This  required  a person  knowledgeable 
of  that  program  to  get  involved  and  greatly  Increased  both  the  risk 
and  the  cost  associated  with  the  simulation  results.  Since  the  thrust 
of  this  total  effort  was  to  provide  an  unchanging  vehicle  dynamics 
package  that  the  engineer  could  use  and  grow  confortable  with,  it  was 
determined  that  the  addition  of  a post  processor  program  to  the  package 
was  required. 

This  program  is  that  post  processor.  A survey  was  made  of  the 
types  of  output  that  users  have  requested  in  the  past.  This  program 
attempts  to  provide  access  to  the  output  forms  desired  for  a wide 
variety  of  signals  from  the  simulation.  Since  it  was  specifically 
developed  to  accept  the  simulation  model  output  included  in  this 
packaoe,  it  was  also  possible  to  write  a user  oriented  interface. 

This  interface  portion  of  the  program  allows  the  user  to  specify 
the  signal  he  desires  in  terms  he  can  understand.  The  signal  is  then 
constructed  as  a linear  combination  of  the  required,  stored  signals 
and  automatically  generated.  This  reduces  the  chance  of  the  user 
Inadvertantly  plotting  the  wrong  signal  and  perhaps  never  knowing  that 
he  did.  As  long  as  the  automatic  signal  generation  capability  is  used, 
the  siqnal  will  be  correctly  labeled  on  all  output  forms  generated. 

The  automatically  generated  signals  that  can  be  selected  Include 
vertical  motion  of  any  point  on  the  sprung  mass  (i.e.  chasis  or  hull), 
pitch  motion  of  sprung  mass  about  Its  CG,  the  vertical  motion  of  any 
wheel  or  the  forcing  function  (terrain  profile)  at  any  wheel,  or  the 
relative  motion  across  any  suspension  system  component.  This  last 
item  corresponds  to  the  spring  deflection  or  shock  absorber  velocity. 

When  the  signal  values  for  this  relative  motion  are  greater  than 
zero,  the  travel  is  in  the  jounce  direction  (I.e.  the  direction  that 
tends  to  compress  the  suspension  system  springs). 
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The  user  also  has  the  opportunity  to  create  his  own  signal  by 
specifying  the  option  OTHER  when  the  signal  position  is  requested. 

If  this  is  specified  the  user  is  asked  for  the  label  he  would  like 
to  have  associated  with  that  signal  and  the  proper  units  for  the  signal. 

He  is  then  permitted  to  construct  a signal  as  any  linear  combination  of 
up  to  (and  including)  three  items  that  are  written  in  the  simulation 
output  file.  This  file  contains  the  time  histories  of  the  accelerations, 
velocities  and  displacements  of  the  sprung  CG  pitch  and  vertical  motions 
as  well  as  for  the  vertical  motion  of  each  wheel  center.  Also  included 
is  the  forcing  function  or  profile  height  at  each  wheel.  Certain  vehicle 
and  terrain  characteristics  are  also  contained  in  the  file  and  are 
used  to  label  the  various  forms  of  output  from  this  program. 

In  order  to  create  his  own  signals  the  user  must  give  the  item 
number  of  the  desired  signal  and  the  scale  factor  by  which  he  desires 
that  signal  to  be  multiplied.  The  item  numbers  correspond  to  the  signal's 
position  within  each  record.  The  arrangement  of  a record  is  shown  in 
Table  1.  The  user  is  responsible  for  making  the  signal  units  agree 
with  the  units  label  he  types  in.  To  accomplish  this  it  should  suffice 
to  remember  that  the  primary  units  of  the  data  in  the  simulation  output 
file  are  feet  and  radians.  Therefore,  if  a velocity  signal  were  desired 
in  inches  per  second,  a scale  factor  of  12  would  have  to  be  supplied 
with  the  appropriate  signal  element  to  convert  from  ft/sec  which  is  the 
standard  units  for  linear  velocity. 

It  is  easy  to  see  (and  easier  to  learn  by  experience)  that  using 
this  option  to  construct  one's  own  signals  requires  a great  deal  of 
care  by  the  user  if  the  integrity  of  the  output  is  to  be  maintained. 

It  is  felt  (and  strongly  hoped)  that  in  most  studies  the  automatic 
signal  generation  capability  will  provide  the  desired  output  signals. 

Once  a signal  has  been  constructed,  either  automatically  or  as 
specified  by  the  user,  several  options  are  available  for  processing 
that  signal.  Each  time  an  option  has  been  completed,  the  user  is 
allowed  to  reselect  from  the  list  of  options.  At  this  point  the 
program  may  he  stopped  or  a new  signal  selected,  The  available 
processing  options  include  listing  certain  signal  statistics,  computing 
the  Power  Spectral  Density  (PSD)  of  the  signal,  plotting  the  signal 
as  a function  of  time,  and  generating  a histogram  from  the  signal. 

Samples  of  several  of  the  processing  option  outputs  are  included  in 
Appendix  I. 
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Item  No. 

Signal  Description 

Units 

1 

CG  Vertical  Acceleration 

ft/sec2 

2 

CG  Vertical  Velocity 

ft/sec 

3 

CG  Vertical  Displacement 

ft 

4 

CG  Pitch  Acceleration 

rad/sec2 

5 

CG  Pitch  Velocity 

rad/ sec 

6 

CG  Pitch  Displacement 

rad 

7 

Wheel  1 Vertical  Acceleration 

ft/sec2 

8 

Wheel  1 Vertical  Velocity 

ft/sec 

9 

Wheel  1 Vertical  Position 

ft 

10 

Profile  Height  Beneath  Wheel  1 

ft 

4N+3 

Wheel 

N Vertical  Acceleration 

ft/sec2 

4N+4 

Wheel 

N Vertical  Velocity 

ft/sec 

4N+5 

Wheel 

N Vertical  Position 

ft 

4N+6 

Profile  Height  Beneath  Wheel  N 

ft 

Table  1 Signal  Position  Within  The  Record  For  The  Binary  Data  File 
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The  full  power  of  this  package  can  only  be  realized  by  using 
a graphics  terminal.  Most  of  the  options,  however,  are  available  in 
some  form  even  if  a graphics  terminal  is  not  used.  Following  comput- 
ation of  either  the  PSD  or  the  histogram  for  a signal  the  user  has 
the  opportunity  to  ask  for  either  a tabular  listing  or  a plot  of  the 
respective  output.  The  statistics  (this  includes  the  number  of  samples, 
minimum,  maximum,  mean  and  rms  levels  of  the  signal)  and  the  vertical 
absorbed  power  do  not  make  use  of  the  graphics.  It  should  be  noted 
that  the  results  of  the  absorbed  power  option  will  have  their  normal 
significance  only  if  the  signal  specified  represents  a vertical 
acceleration  in  ft/sec^.  The  signal  plotting  option  of  course  cannot 
be  employed  unless  the  program  is  being  run  from  a graphics  terminal. 

A local  file  name  PRINT  will  be  generated  whenever  this  program 
is  run.  This  file  will  contain  at  least  limited  results  for  each 
option  that  is  run.  This  will  always  include  certain  simulation  and 
signal  identification  followed  by  a limited  set  of  signals  statistics 
(the  rms  will  not  be  Included  unless  the  statistics,  PSD,  or  absorbed 
power  options  are  run).  If  the  PSD  option  is  run  and  tabular  output 
is  specified,  a listing  of  signal  power  and  signal  amplitude  as  a function 
of  frequency  will  bo  written  in  this  PRINT  file.  In  this  case  rio 
similar  listing  will  be  presented  at  the  user's  terminal.  This  is 
the  only  situation  In  running  this  program  where  this  is  done.  The 
rationale  for  this  is  tnat  these  listings  may  be  as  long  as  ZOO  lines 
and  the  time  sharing  terminals  are  very  slow.  The  contents  of  the 
PRINT  file  may  be  obtained  in  either  of  two  ways.  Once  the  user  has 
stopped  this  program  he  may  still  list  the  PRINT  file  at  his  terminal 
with  the  following  two  statements: 

REWIND,  PRINT 

COPY,  PRINT,  OUTPUT 

■f  instead  it  is  desired  to  print  the  file  on  a remote  batch 
terminal  the  following  two  statements  will  suffice: 

REWIND,  PRINT 

HATCH,  PRINT,  PRINT, B5,  NAME 


?Z 


In  the  BATCH  command,  the  first  occurence  of  PRINT  refers  to  the 
file  name.  The  parameter  B5  Is  the  terminal  ID  at  which  the  output 
Is  to  be  printed  (B5  is  currently  the  ID  for  the  batch  terminal  at 
TARADCOM).  The  final  parameter  NAME  is  some  arbitrary  name  that  will 
be  associated  with  the  printout. 

The  running  of  this  post  processor  program  (l.e.  BANALYZE)  can  be 
Initiated  with  the  following  two  statements: 

ATTACH , BAN,  BANALYZE,  ID-DRDTAZSA 

BAN,  BINARY 

The  parameter  BINARY  In  the  second  statement  above,  should  be  re- 
placed with  the  local  file  name  corresponding  to  the  binary  data  file 
which  Is  to  be  analyzed.  This  file  should  have  been  created  by  the 
simulation  package  described  within  this  report. 

It  is  necessary,  of  course,  to  run  this  program  from  a graphics 
terminal  {a  Tektronix  4014)  If  any  plotting  Is  to  be  done.  If  the 
plotting  options  are  not  used,  any  Interactive  terminal  will  be 
adequate. 
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One  type  of  vehicle  dynamic  analysis  often  requested  is  to  find 
the  driver's  ride  limiting  velocity  across  various  cross  country 
terrains.  This  analysis  can  be  performed  using  the  program  BVRIDE. 

The  problem  also  arises  to  determine  the  driver's  ride  limiting  speed 
over  discrete  semicircular  obstacles  of  various  heights.  The  ride 
limiting  criteria  for  the  driver  over  discrete  cbstacles  is  usually 
considered  to  be  2 . 5q  peak  vertical  acceleration,  whereas  for  cross 
country  driving  an  overage  vertical  absorbed  power  limit  of  six  watts 
is  used.  The  BVRIDE  program  prints  the  peak  acceleration  and  average 
absorbed  power  for  each  simulation  and  could  be  used  for  the  ride 
performance  evaluations  over  discrete  obstacles. 

However,  since  several  simulations  are  run  for  each  obstacle  height 
and  each  one  is  quite  short,  a modified  version  of  BVRIDE  is  offered 
for  simulating  the  vehicle  over  semicircular  obstacles.  The  advantage 
offered  by  this  modified  program  is  that  up  to  ten  vehicle  velocities 
can  be  run  over  a particular  user  specified  obstacle  with  a single 
run  command. 

The  output  is  a table  of  vehicle  velocity  versus  the  simulated 
driver's  peak  vertical  acceleration  at  each  velocity.  Rather  than 
requiring  a single  velocity  the  user  must  supply  an  initial  velocity, 
a delta  velocity,  and  a final  velocity.  Runs  will  be  made  at  increasing 
velocities  until  one  of  three  things  occur.  The  three  possible  con- 
ditions that  will  cause  the  run  to  terminate  are:  the  maximum  velocity 
requested  has  been  exceeded;  ten  different  velocities  have  been  simulated; 
or  a run  has  been  made  in  which  the  driver  experienced  at  least  2 - 5g 
vertical  acceleration. 

This  program  acepts  the  same  vehicle  data  file  as  does  the  program 
f’-VIDE  and  the  time  limit  may  have  to  be  extended  for  this  program, 
as  shown  in  the  RVRIDF  description.  No  terrain  data  file  is  require, 
however,  since  the  appropriate  obstacle  is  generated  within  this  program. 

A binary  file  can  be  generated  by  this  program  which  can  he  analyzed 
with  the  post  processor  program  BANALYZE.  This  binary  file  will  not 
normally  be  created  but  when  it  is  it  will  contain  the  results  from 
a single  simulation  (i.e,  at  only  one  vehicle  velocity).  If  a binary 
file  is  desired,  the  user  can  cause  its  generation  by  entering  a zero 
for  the  maximum  vehicle  velocity.  If  this  is  done  a single  simulation 
will  be  performed  at  the  minimum  velocity  entered,  the  results  will  be 
printed  and  the  binary  output  file  will  be  written. 
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time  sharing  terminal  with  the 


This  program  can  be  run  from  any 
following  two  statements: 

ATTACH,  ESWHEELS,  I D*DRDTAZ$A 

B WHEELS, VEH, BINARY 

The  parameters  VEH  and  BINARY  have  the  same  interpretation  here  as 
they  do  in  the  running  of  BVRIDE. 


Input  / 

Vehicle  & / 

\ Terrain  / 

\ Data  j subroutine 


-T  UFA 


Acceleration 
this  wheel 


c.g.  vertical  and 
pitch  accelerations 


LINEAR 


YCALC 


function 

calculate  wheel  path 


initialize  integration  routines 
initialize  vehicle 

get  wheel  position 


integrate,  wheel  equations 


tandom  equations 
(if  present) 


integrate  c.g.  equations 

driver's  acceleration 


update  integration  package 


driver's  absorbed  power 


VRTDF.  Flow  Diagram  by  Subroutine 
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APPENDIX  I 

Sample  Vehicle  Analysis 
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Engineering  Data  Required  for  Test  Jeep 


Sprung  Weight 2710  lbs 

Unsprung  Weight 410  lbs 

Sprung  Pitch  Moment  of  Inertia  about 

the  CG  719  slug-ft2 

Wheel  No.  1 Forward  of  CG 44  in 

Wheel  No.  2 Aft  of  CG 40  In 

Driver's  Position  Forward  of  CG 1 In 

Spring  Rate  of  Tires  675  lbs/1  n 

Damping  Rate  of  Tires  0.75  lbs-sec/ln 

The  spring  and  damping  curves  representing  the  vehicles  suspension 
system  are  shown  In  the  following  two  pages. 

This  test  vehicle  is  similar  to  the  data  for  the  M151A1  1/4  ton  truck. 
Certain  minor  changes  were  made,  however,  In  the  interest  of  simplifying 
the  data  Input  example.  If  precise  data  on  the  M151A1  Is  desired,  It 
should  be  obtained  from  another  source. 
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CURVE  FOR  THE  WHEEL  SPRINGS 
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ICAL  DISPLACEMENT 


fUSlftl  i/4  TOW  TRUCK 

10.#  WPH  THREE  INCH  SPACED  KJRP  COURSE 


FREQUENCY — H2 


APPENDIX  II 


Vehicle  Data  File  Contents 


CARD  1 


FORMAT  (6A10) 


The  first  60  columns  of  this  card  arc  used  for  the  vehicle  name  and 
description.  This  information  will  serve  not  only  to  identify  the  file 
but  will  also  appear  on  all  outputs  from  this  package. 


CARD  2 

FORMAT  (5F10.0,  412) 

WS 

Sprung  weight  of  vehicle  (lbs) 

wu 

Unsprung  weight  of  vehicle  (lbs) 

J 

Sprung  pitch  moment  of  inertia  of  vehicle  (slug-ft^) 

XC(1) 

Horizontal  distance  of  selected  location  from  CG  (ft). 

This  must  be  negative  if  location  is  behind  the  CG. 

YC(1) 

Vertical  distance  of  selected  location  from  CG  (ft). 

This  must  be  negative  if  location  is  CG.  Note:  this 
location  is  referred  to  as  the  driver's  position  on 
output . 

N 

Number  of  wheels  on  one  side  of  vehicle. 

NCG 

Number  of  wheels  in  front  of  CG  on  one  side  of  vehicle. 

NT 

Number  of  tandems  on  one  side  of  vehicle. 

NTCG 

Number  of  tandems  in  front  of  CG  on  one  side  of  vehicle. 
The  tandem  is  considered  to  be  in  front  of  the  CG,  if, 
and  only  if,  the  pivot  point  of  the  tandem  is  in  front 
of  the  CG. 

CMd)  3 

FORMAT  (7F10.0) 

J,  ( J l 

For  I = 1 to  N,  the  number  of  wheels  per  side.  Horizontal 
distance  from  CG  to  Ith  wheel.  Enter  absolute  value  of 
number  (ft).  If  there  are  more  than  seven  wheels,  on  a 
side,  two  cards  will  be  needed. 

CARD  4 

FORMAT  (6011) 

TNDARY 

(I) 

Enter  a 0 if  the  wheel  is  not  in  tandem,  a 1 if  it  is  the 
front  wheel  in  a tandem,  and  a 2 if  it  is  the  second  wheel 
in  a tandem. 

NOSSEG 

(l) 

Number  of  segments  in  vehicle's  spring  curves,  (maximum 
is  6) . 

nodshg 

(I) 

Number  of  segments  in  vehicle's  damping  curves,  (maximum 
is  6) . 

3° 


DMPTYP 

NDRSEG 

NORDEG 

NOTE: 

CARD  5 
VELOPT 

POOPT 

SPOFT 

DAMDl’T 

SRUBOP 

DRUBOP 

KTYPli 

NLTrii: 


(I)  Enter  a 0 to  indicate  no  damper;  a 1 for  a one-way  damper 

(i.e. , damping  on  the  compression  stroke  only);  and  a 2 
for  a two-way  damper. 

(I)  Number  of  segments  in  tire's  spring  curve,  (maximum  is  6). 

(I)  Number  of  segments  in  tire's  damping  curve,  (maximum  is 

6). 

I goes  from  1 up  through  N,  the  number  of  wheels  on  one  side.  The 
N values  of  one  variable  are  entered  before  any  of  the  next  variable. 

FORMAT  (714) 

Inter  a 1 if  a sprocket  is  to  be  included  in  the  simulation, 
otherwise  a 0.  If  a one  is  entered,  CARD  12  must  be  included. 

Enter  a 0 if  the  unsprung  weight  is  to  be  divided  evenly 
among  the  wheels  or  a 1 if  the  weight  of  each  wheel  is  to 
be  entered.  If  a 1 is  entered,  then  CARD  7 must  be 
present. 

A 1 signifies  all  vehicle's  spring  curves  are  the  same  (thus 
only  the  first  need  be  in  the  data  deck  providing  NOSSEG 
(I)  = 0 on  CARD  5 for  II).  A value  greater  than  1 indicates 
the  number  of  cures  to  be  read  in.  The  Cards  8 thru  10 
are  needed  for  each  wheel  that  has  NOSSEG  (I)  0.  The 

spring  curve  for  any  wheel  that  is  not  read  in  will  be  set 
the  same  as  the  curve  of  the  previous  wheel. 

A 1 signifies  all  vehicle  damping  curves  are  the  same. 

A 2 indicates  that  all  the  wheels  do  not  have  the  same 
damping.  If  only  certain  wheels  have  damping,  a 2 must 
be  entered. 

A 1 indicates  all  tires  have  the  same  spring  curve;  a 2 
means  they  don't  all  liave  the  same  spring  curve. 

A 1 indicates  that  all  tires  have  the  same  damping  curve; 
a 2 means  they  don't  all  have  the  same  damping  curve. 

Enter  a 1 if  this  vehicle  has  a fifth  wheel  hitch,  otherwise 
enter  a 0.  If  a 1 isentcred,  CARD  13  is  required. 

At  present,  at  least  one  segment  must  be  used  for  the  tire,  spring, 
and  damping  curves. 
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CARP  6 


FORMAT  (7  F10.0) 


This  card  is  optional • It  is  needed  only  when  the  value  of  POOPT  on 
CARP  5 is  greater  than  zero.  The  sum  of  the  wheel  weights  should  he  the 
unsprung  weight  for  one  side  of  the  vehicle. 

MW  (I)  For  T**l  to  N.  The  weight  of  the  Ith  wheel  (lbs). 


CARP  7 FORMAT  (3  F10.0) 


One  of  these  cards  must  be  supplied  for  each  tandem  I for  I«*l,  NT. 


TNPJ  (I) 


2 

The  pitch  moment  of  inertia  of  the  tandem  (slug-ft  ) 
about  its  pivot  point. 


l.WB  (I) 


Distance  from  wheel  center  to  tandem  pivot  (feet). 


LK  (T) 


Horizontal  distance  from  vehicle's  sprung  CC  to  tandem 
pivot  (feet) . 


CARPS  8 - 10  FORMAT  (6  F10.0) 

This  set  of  data  cards  will  come  in  groups  of  three.  These  cards  will  be 
used  to  describe  the  curves  for  the  various  spring  and  damper  elements 
of  the  vehicle  model. 


The  values  of  NOSSF.C  (I),  NOPSRO  (I),  NORSRO  (I),  AMP  MORPF.O  (I)  were 
obtained  from  CARP  4 for  each  wheel  T.  These  values  will  then  be  checked 
in  that  order  for  each  wheel . Each  time  a positive  value  is  found,  call 
that  value  Ml , three  cards  will  he  read  in.  The  first  card  will  contain 
the  El  slopes  of  the  linear  segments  of  the  curve.  The  second  card  will 
contain  the  Y-intcrcepts  of  each  of  these  N1  segments  of  the  curve.  The 
third  card  will  then  contain  the  N1  horizontal  values  at  which  the  corres- 
ponding segments  end.  The  values  on  these  three  cards  should  be  in  LBS/FT, 
respectively  (or  I.BS-SF.C/FT,  LBS,  and  FT/SEC  for  damping  curves). 


Whenever  n 0 is  found  for  the  number  of  segments  of  a curve,  no  data 
cards  are  rend  In.  This  process  is  continued  for  I-l  up  to  and  including 
T=N  where  N is  the  number  of  wheels  on  one  side  of  the  vehicle. 


Tt  should  he  noted  that  if  the  vehicle  spring  curves,  vehicle  damping 
curves,  tire  spring  curves,  or  tire  damping  curves  are  identical  for 
all  wheels,  then  only  the  curve  for  the  first  wheel  needs  to  be  read  in. 
In  order  to  do  this,  however,  the  number  of  segments  for  that  curve  must 


dl 


be  equal  to  zero  for  all  other  wheels  on  CARD  5.  Also,  the  appropriate 
option  must  be  set  equal  to  one  on  CARD  5 (SPOPT,  DAMOPT,  SRUBOP  or 
DRIJBOP) . 

It  should  be  noted  also  that  all  curves  needed  for  one  wheel  must 
be  supplied  before  any  curves  for  the  next  wheel. 

/ 

This  program  determines  the  average  static  load  required  at  each 
wheel  to  support  the  vehicle.  The  actual  static  position  Is  then 
assumed  to  be  on  the  same  linear  segment  as  the  average  static  load 
is.  In  other  words,  once  the  proper  segment  is  determined,  this 
segment  is  assumed  to  be  extended  over  the  necessary  range  of  deflection 
in  order  to  carry  out  the  initialization  of  the  vehicle.  Once  the 
initial  conditions  are  determined  the  actual  spring  curve  representations 
are  used  in  the  simulation. 


CARD  11  FORMAT  (F10.0) 

DIAM  The  effective  wheel  diameter  (i.e.,  twice  the  rolling 

radius)  in  inches. 


CARD  12  FORMAT  (3F10.0) 

This  card  is  needed  only  if  the  value  of  VELOPT  on  CARD  5 is  equal 
to  1.  This  card  provides  the  necessary  data  on  the  sprocket  if  one 
is  to  be  Included  in  the  simulation. 

SPRCL  The  vertical  sprocket  clearance  at  static  (ft). 

SPROC  The  horizontal  distance  of  the  sprocket  from  the 
CG  (ft). 

SPRK  The  spring  rate  of  the  sprocket  (lbs/ft). 
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CARD  13 


FORMAT  (3  FIO.O) 


This  card  is  needed  if  the  value  of  KTYPE  on  CARD  5 is  equal  to  1. 

This  card  provides  the  necessary  data  for  the  hitch  on  a truck  tractor 
when  modeling  tractor-trailer  combinations. 

CGHT  Horizontal  distance  from  vehicle  CG  to  fifth  wheel 

hitch  (ft). 

SRHT  Spring  rate  of  fifth  wheel  hitch  (Ibs/ft). 

DRHT  Damping  rate  of  fifth  wheel  hitch  (Ibs/ft/sec) . 
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APPENDIX  III 

Terrain  Data  File  Format 
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CARD  1 
NTER 

D2SINC 

TERIDT 

CARD  2 
ATAB 


FORMAT  (I5,F10.0,5A10) 

The  number  of  terrain  data  points  that  are  to  be  read 
in. 

The  hori zonal  spacing  of  the  terrain  data  points  (in  feet). 

A terrain  label  of  up  to  50  characters  in  length  may  be 
entered  here.  This  label  v/ill  appear  on  all  printouts 
using  this  terrain. 


thru  LAST  'CARD 


FORMAT  (6F10.0) 

The  terrain  elevations  in  feet  at  the 
specified  hori zonal  spacing. 
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appendix  IV 
Standard  Terrain  Files 


4f» 


' 

Terrain  catalog 

FOR  DROTAZSA 

FILE  NAME 

RM5TIN, ) 

LENGTH  IPTS/i 

DESCRIPTION  (DAtE) 

APG9 

1.04  RMS 

301  FT.  (301) 

APG  COURSE  NINE 

aPGTI 

APG12 

1732  ~Rm§ 
1.70  RMS 

401  FIT  (401) 

301  FT.  (301) 

APG  COURSE  ELEVEN 

APG  COURSE  TWELVE 

~msn 

"APirer 

~AFG3r 

APG3TI 

APG37 

APGIV 


3~0T~FTY~“<  30 1) 

^5rrrr*<25i) 


T.  1 5TFH5” 


jp~s  course  twent y-Ffi net 

(11  NOV  76) 

' APG~COOR3FTHTfm 

(10  NOV  76) 


“3T£nsn?R5“ 
‘0 .66“ RMS" 
T.< 


TroTTTrrroTr 
T01~FTV~OTrr 
1469) 


"A  P G “COU  r5F~T  RIRTY^'CWE — 
(16  NOV  76) 

APG"COGRSrTffIftTY-SEWN 
<13  NOV  76) 


— fthdti sTomws nnrrT.  TSTooi rr;"flooo“TRTrm 

FTHDT3  0.60  RMS  401  FT.  (401)  FT.  H000  TRAIL  T3 

— FTHDcciA — tn&trpm 4W"rrr~r40O7 ft.~r  qoo~tratc— ccta 

FTHOSR1  0.40  RMS  401  FT.  (401)  FT.  HQQO  SECONDARY  ROAD 


TERR4 1.94  RMS  372  FT,  (372)  FT.  KNOX  ME01UH  TERRAIN 

PERRY 4.12  RMS  442  FT.  (44?)  PERRYMAN  TERRAIN  SIX 

LEVELRD  0.00  RMS  132  FT.  (132)  FLAT  LEVEL  ROAD  (30  AUG  76) 


BUMP6  X.XX  RMS  95  FT.  (190)  ONE  6 INCH  TRAPEZOIDAL  BUMP 

— buhps3 xtxxhrms iwm  rrwi HWErnwnrTFAi^^  bumps 

SPACED  31  FT,  APART. 

• - R0UND3  -Jr.XK“RMS  I BTXTTT  T22TH  OKiF~TTNCH  ^EMITTRCULIR  MUMP  ' 

ROUND6  X.XX  RMS  11  FT.  (132)  ONE  6 INCH  SEMICIRCULAR  BUMP 

ROUNDS X.KX-RMS B4'  FT._T370) ONE"fl  INCF'T5CTTClFCt)PF_Bl3MP~ 

ROUN012  X.XX  RMS  70.5  FT.  (370)  ONE  12  INCH  SEMICIRCULAR  BUMP 


FTNX01  .75  RMS  501  FT.  (501)  FT.  KNOX  COURSE  Nl  RT,  TRACK 

fFnHIfrTB1 

FTNXQ2  .84  RMS  501  FT.  (501)  FT.  KNOX  COURSE  61  LEFT  TRACK 

— fgg  “"HURT  1 TB ) ‘ " — — 

FTNX03  1.10  RMS  501  FT.  (501)  FT.  KNOX  COURSE  6*  RT.  TRACK 

r?4~MAR-7a) 


FTNX04 

1.24  RMS 

501 

FT. 

(501) 

FT.  KNOX  COURSE 

62 

LEFT 

TRACi* 

FTNX05 

1.80  RMS 

501 

FT. 

(501), 

■ -nfl~HAir~7A) 

FT.  KNOX  COURSE 

63 

LEET 

TRUCK 

FTNX06 

1.95  RMS 

501 

FT. 

(501) 

(23^WAR~78r_ 

FT.  KNOX  COURSE 

63 

RT. 

trac» 

FTNX07 

2.84  RMS 

401 

FJ, 

(401) 

rzTin&rTKi 

FT.  KNOX  COURSE 

64 

_LFFT 

TRACK 

(13-MART1T 
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terrain  catalog  for  ordtazsa  pagf  i 


FILE  NAME  RMS (IN, ) LENGTH  (RTS.)  DESCRIPTION  (DATE) 


FTNX08 

3.61 

RMS 

401 

FT. 

(401) 

FT. 

KNOX  COURSE 

#4 

RT.  TRACK 

FTNX09 

5.22 

RMS' 

401 

FT. 

(401) 

(IB 

FT. 

MAR  78) 

KNOX  COURSE 

#5 

RT,  TRACK 

FTNXIO 

5.27 

RMS 

401 

FT. 

(401) 

(16 

FT. 

MAR  78) 

KNOX  COURSE 

MS 

LEFT  TRACK 

FTKN1 

.82 

RMS 

501 

FT. 

(501) 

(18 

FT. 

MAR  76) 

KNOX  COURSE 

#1 

RT.  TRACK 

FTKN2 

.83 

RMS 

501 

FT. 

(501) 

( 5 
FT. 

MAY  78) 

KNOX  COURSE 

Ml 

LEFT  TRACK 

FTKN3 

1.03 

RMS 

501 

FT. 

<5015 

(27 

FT. 

MAY  78) 

KNOX  COURSE 

M2 

RT.  TRACK 

FTKN4 

1.04 

RMS 

501 

FT. 

(501) 

( 5 
FT. 

MAY  78) 

KNOX  COURSE 

M2 

LEFT  TRACK 

FTKNS 

2.07 

RMS 

501 

FT, 

(SOI) 

<26 

FT. 

MAY  78) 

KNOX  COURSE 

M3 

RT.  TRACK 

FTKN6 

1.85 

RMS 

501 

FT. 

(501) 

( 1 
FT. 

MAY  78) 

KNOX  COURSE 

M3 

LEFT  TRACK 

FTKN7 

3.91 

RMS 

401 

FT, 

(401) 

( 2 
FT. 

MAY  78) 

KNOX  COURSE 

M4 

RT.  TRACK 

FTKN8 

3. 57 

RMS 

401 

FT. 

(401) 

( 3 

FT. 

MAY  78) 

KNOX  COURSE 

M4 

LEFT  TRACK 

FTKN9 

2.31 

RMS 

401 

FT. 

(401) 

( 5 MAY  78) 

FT.  KNOX  COURSE 

M5 

RT.  TRACK 

FTKNIO 

2.40 

RMS 

401 

FT, 

(401) 

( 3 
FT. 

MAY  78) 

KNOX  COURSE 

MS 

LEFT  TRACK 

( 1 MAY  18) 

FORT2 

.87 

RMiT 

^01 

nr 

(501) 

FT. 

KNOX  COURSE 

Ml 

YEFT  TRACK 

.99 

(14 

JUL  78) 

FORTS 

RM5~~ 

501  FT, 

(501) 

fTT 

<14 

KNOX  COURSE 
JUL  78) 

M2 

rYY1rack’~ 

FORT4 

1.05 

RMS 

501 

FT. 

(501) 

FT.  KNOX COURSE 

M2 

left  Yrack 

FORTS 

<14 

JUL  76) 

176^ 

“R  M5 

"Strr 

TT7 

rsoT) 

FT  . 

KNOX  COURSE  M 3 RT.  TRACK 

501 

<14 

JUL  78) 

m3 

F0RT6 

1.88  RMS 

FTV  (501) 

FT.  KNOX  COURSE 
<14  JUL  78) 

LeFT  track 

F0RT7 

1 .T3KM5 

401  TT^C401Y 

FT.  KNOX  COURSr  #41?T.  I'RACK 

(14 

JUL  78) 

FOR  1 8 

3 * F 4 

"RMS 

~i nrr 

FT, 

"iron 

FT. 

knox  CouRSe 

K4~ 

IETTWCK  ’ 

(14 

JUl.  78) 

FORT  9 

2.28 

__  g > 4T 

RMS 

4 01 

FT.' 

(401) 

' FT. 

<14 

KNOx  COURSE 
JUL  78) 

#5  rt . Track 

FORT 10 

~RM5 

401 

PT  9 

14011 

'FT. 

KNOx  'COURSE 

M5 

LEFTThaCK 

( 1 A JUL  78) 


48 


DISTRIBUTION  LIST 


Commander 

US  Army  Tank-Automotive  Research  A Development 
Command 

ATTN:  Applied  Research  Function,  DRDTA-ZSA 

Survival  Technoloqy  Function,  DRDTA-ZSS 
Technical  Library,  DRDTA-UL 
Warren,  MI  48090 

Commander 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 


SECURITY  CL  ASM  VACATION  or  THIS  PAOE  !»*•«  D *•  Bnlorod) 


|T  rEpor?  number 

1241 3^ 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  OOVT  ACCESSION  NO 


1.  RECIPIENT'S  CATALOG  HUMBER 


4.  TITLE  (mtd  SMIli) 


».  TYPE  Of  REPORT  * PERIOD  COVERED 


Interactive  Vehicle  Dynamics  and  Ride 
Evaluation  Package 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!*; 

Francis  B.  Hoogterp 


I.  CONTRACT  OR  GRANT  NUMBER!*; 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Tk-Autmv  Rsch  & Dev  Cornnand 
DRDTA-ZSA  Warren,  MI  48090 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  • WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

US  Am\y  Tk-Autmv  Rsch  & Dev  Command 
DRDTA-UL 

Warren,  MI  48090 

14.  MONITORING  AGENCY~NAME  4 ADDRESS!/!  dltloro nl  from  Cetilrolllnd  Olll oaj 


12.  REPORT  DATE 


11.  NUMBER  OF  PAGES 

50 


IS.  SECURITY  CLASS,  (ol  I hi*  r*portJ 

UNCLASSIFIED 


1»-.  DECL  ASM  FFCATtONA>©WN  GRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ol  thlt  Ropori) 


Approved  for  public  release,  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  tho  abtlraul  onlorod  In  B look  JO,  II  dllloront  hom  Ropori) 


It.  SUPPLEMENTARY  NOTES 


IS.  KEV  WORDS  (Conti nut  on  rorotoo  tldo  II  nmooooory  and  Idontlly  by  blook  numbor) 

Vehicle  Dynamics 
Ride  Evaluation 


20.  ABSTRACT  (Continue  oft  r*var««  «Jt/«  U ttmcmmm^ry  mtd  identity  by  bJoclr  numb+f) 

J^Thls  report  describes  an  Integrated  package  for  evaluating  the  dynamic  effects 
^ of  driving  a vehicle  across  any  described  section  of  terrain.  The  model  in- 
corporated includes  only  the  pitch  and  roll  motions  of  the  chasls  (or  hull). 

The  package  Is  composed  of  a four  digital  computer  programs.  Each  program  is 
designed  to  be  run  Interactively  and  does  not  require  that  the  user  have  any 
programming  experience.  The  programs  provide  the  capability  to  create  vehicle 
data  files,  perform  simulations,  and  specify  the  form  and  types  of  the  output 
desired.  A general  description  of  the  model  structure,  its  capabilities  and  its 
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